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Abstract 

One still cannot conclusively assert that the universe underwent a 
neutral phase, despite the new COBE FIRAS hmit y < 2.5 x 10^^ on 
Compton y-distortions of the cosmic microwave background. Although 
scenarios where the very early {z ~ 1000) ionization is thermal (caused 
by IGM temperatures exceeding lO'^K) are clearly ruled out, there is 
a significant loophole for cosmologies with typical CDM parameters 
if the dominant ionization mechanism is photoionization. If the ion- 
izing radiation has a typical quasar spectrum, then the y-constraint 
implies roughly h^^'^Qigm^Q^''^^ < 0.06 for fully ionized models. This 
means that BDM models with fio ~ 0.15 and reionization at z « 1000 
are strongly constrained even in this very conservative case, and can 
survive the y test only if most of the baryons form BDM around the 
reionization epoch. 
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1 Introduction 



Recombination of the primeval plasma is commonly assumed but was by 
no means inevitable. Theories exist that predict early reionization are as 
diverse as those invoking primordial seed fluctuations that underwent early 
collapse and generated sources of ionizing radiation, and models involving 
decaying or annihilating particles. The former class includes cosmic strings 
and textures, as well as primordial isocurvaturc baryon fluctuations. The 
latter category includes baryon symmetric cosmologies as well as decaying 
exotic particles or neutrinos. 

The Compton ^/-distortion of the cosmic microwave background (CBR) 
provides a unique constraint on the epoch of reionization. In view of the 
extremely sensitive recent FIRAS limit of y < 2.5 x 10~^, we have reinvesti- 
gated constraints on the early ionization history of the intergalactic medium 
(IGM), and have chosen to focus on what we regard as the most important 
of the non-standard recombination history models, namely the primordial 
isocurvature baryon scenario involving a universe dominated by baryonic 
dark matter (BDM), as advocated by Peebles (1987); Gnedin & Ostriker 
(1992) (hereafter "GO"); Cen, Ostriker & Peebles (1993) and others. This 
class of models takes the simplest matter content for the universe, namely 
baryons, to constitute dark matter in an amount that is directly observed 
and is even within the bounds of primordial nucleosynthesis, if interpreted 
liberally, and can reconstruct essentially all of the observed phenomena that 
constrain large-scale structure. The BDM model is a non-starter unless the 
IGM underwent very early reionization, in order to avoid producing exces- 
sive CBR fluctuations on degree scales. Fortunately, early nonlinearity is 
inevitable with BDM initial conditions, 5p/ p oc M~^/^^, corresponding to 
a power-spectrum {5\) oc k~^^'^ for the obscrvationally preferred choice of 
spectral index (Cen, Ostriker & Peebles 1993). 

Is it possible that the IGM has been highly ionized since close to the 
standard recombination epoch at z ~ 1100? Perhaps the most carefully 
studied BDM scenario in which this happens is that by GO. In their sce- 
nario, Oq = r^bo ~ 0.15. Shortly after recombination, a large fraction of 
the mass condenses into faint stars or massive black holes, releasing energy 
that reionizes the universe and heats it to T > 10, OOOK by z = 800, so 
Compton scattering off of hot electrons causes strong spectral distortions 
in the cosmic microwave background. The models in GO give a Compton 
y-parameter between 0.96 x 10""^ and 3.1 x 10~^, and are thus all ruled 
out by the most recent observational constraint from the COBE FIRAS 
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experiment, y < 2.5 x 10~^ (Mather et al. 1994). 

There are essentiahy four mechanisms that can heat the IGM sufficiently 
to produce Compton y-distortions: 

• Photoionization heating from UV photons (Shapiro &: Giroux 1987; 
Donahue & Shull 1991) 

• Compton heating from UV photons 

• Mechanical heating from supernova-driven winds (Schwartz et al. 1975; 

Ikcuchi 1981; Ostriker k Cowic 1981) 

• Cosmic ray heating (Ginzburg &; Ozernoi 1965) 

The second effect tends to drive the IGM temperature towards two-thirds 
of the temperature of the ionizing radiation, whereas the first effect tends to 
drive the temperature towards a lower value T* that will be defined below. 
The third and fourth effect can produce much higher temperatures, often 
in the millions of degrees. The higher the temperature, the greater the 
y-distortion. 

In the GO models, the second effect dominates, which is why they fail 
so badly. In this paper, we wish to place limits that are virtually impossible 
to evade. Thus we will use the most cautions assumptions possible, and 
assume that the latter three heating mechanisms are negligible. 



2 The Compton ^-Parameter 

Thomson scattering between CBR photons and hot electrons affects the 
spectrum of the CBR. It has long been known that hot ionized IGM causes 
spectral distortions to the CBR, known as the Sunyaev-Zel'dovich effect. A 
useful measure of this distortion is the Comptonization y-parameter (Kom- 
paneets 1957; Zel'dovich & Sunyaev 1969; Stebbins & Silk 1986; Bartlett & 
Stebbins 1991) 



y 



2 



UeCTtc dt = y* 



jl + z) 



ATi{z)x{z)dz, (1) 



where 



y = 



XHe 

Ax 



Y 



k X 10^K\ / 3Hoi},igmCrtc 



rUeC^ 



SirGrrin 



9.58 X W^ha 



igm- 



2 



Here Tg is the electron temperature, is the CBR temperature, AT4 = 
(Te — T^)/10'^K, ^igm is the fraction of critical density in intergalactic 
medium, and x{z) is the fraction of the hydrogen that is ionized at red- 
shift z. Note that we may have Vlig^ ^ ^-e. all baryons may not be 
in diffuse form. The integral is to be taken from the reionization epoch 
to today. In estimating the electron density ng, we have taken the mass 
fraction of helium to be F « 24% and assumed XHe ~ a;, i.e. that helium 
never becomes doubly ionized and that the fraction that is singly ionized 
equals the fraction of hydrogen that is ionized. The latter is a very crude 
approximation, but makes a difference of only 6%. 

Let us estimate this integral by making the approximation that the IGM 
is cold and neutral until a redshift Zjon, at which it suddenly becomes ionized, 
and after which it remains completely ionized with a constant temperature 
Tg. Then for Zic^n 3> 1 and Tg » Zicm x 2.7K we obtain 

y «i 6.4 X lQ-^h9.ig^%^'^n z-i^, 

where T4 = Tg/IO^K. Substituting the most recent observational constraint 
from the COBE FIRAS experiment, y < 2.5 x 10"^ (Mather et al. 1994), 
into this expression yields 

..„.< 554T-/^fiy^(^)"''. (2) 

Thus the only way to have Zion as high as 1100 is to have temperatures 
considerably below lO^K. In the following section, we will see to what 
extent this is possible. 

3 IGM Evolution in the Strong UV Flux Limit 

In this section, we will calculate the thermal evolution of IGM for which 

• the IGM remains almost completely ionized at all times, 

• the Compton y-distortion is minimized given this constraint. 

3.1 The ionization fraction 

In a homogeneous IGM at temperature T exposed to a density of ( UV 
photons of energy hv > 13.6 eV per proton, the ionization fraction x evolves 
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as follows: 



dx 1 + z 



d{-z) VI + 



\pi{l - x)+ Xcix{l - x) - XrecX , (3) 



where Hq ^(1 + z) ^ times the rates per baryon for photoionization, colli- 
sional ionization and recombination are given by 

Xpi « 1.04 X 10^2 [/lf],g„ai8] C, 

Xci « 2.03 X 104/iJ7i5^r2/2g-i5.8/r4, 



Xrec ^ 0.717/ir2j(^fn,2^4 



1.808 - 0.51nr4 + O.ISTT^''^ 



and T4 = Te/IO^K. Here ais is the spectrally-averaged photoionization 
cross section in units of 10~^®cm2. The differential cross section is given by 
(Osterbrock 1974) 

dais f° ifz^<13.6eV, 

1 6.30^^j7^-:32w7Y ifi^>13.6eV, 



1/4(1-6-2'^/^) 



where 



13.6 eV 



1. 



The recombination rate is the total to all hydrogenic levels (Seaton 1959; 
Spitzer 1968). Recombinations directly to the ground state should be in- 
cluded here, since as will become evident below, the resulting UV photons 
are outnumbered by the UV photons that keep the IGM photoionized in 
the first place, and thus can be neglected when determining the equilibrium 
temperature. 

At high redshifts, the ionization and recombination rates greatly exceed 
the expansion rate of the universe, and the ionization level quickly adjusts to 
a quasi-static equilibrium value for which the expression in square brackets 
in equation (^) vanishes. In the absence of photoionization, an ionization 
fraction x close to unity requires Tg > 15, OOOK. Substituting this into 
equation (|2|) gives consistency with Zion > 1000 only if hQigm < 0.008, 
a value clearly inconsistent with the standard nucleosynthesis constraints 
(Smith et al. 1993). Thus any reheating scenario that relies on collisional 
ionization to keep the IGM ionized at all times may be considered ruled out 
by the COBE FIRAS data. 
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However, this does not rule out all ionized universe scenarios, since pho- 
toionization can achieve the same ionization history while causing a much 
smaller y-distortion. The lowest temperatures (and hence the smallest y- 
distortions) compatible with high ionization will be obtained when the ion- 
izing flux is so strong that Apj ^> Xd- In this limit, to a good approximation, 
equation (^) can be replaced by the following simple model for the IGM: 

• It is completely ionized (x = 1). 

• When a neutral hydrogen atom is formed through recombination, it is 
instantly photoionized again. 

Thus the only unknown parameter is the IGM temperature Tg, which de- 
termines the recombination rate, which in turn equals the photoionization 
rate and thus determines the rate of heating. 

3.2 The spectral parameter T* 

The net effect of a recombination and subsequent photoionization is to re- 
move the kinetic energy ^kT from the plasma and replace it with the kinetic 
energy f/cT*, where T* is defined by ^kT* = {E^v) - 13.6 eV and is 
the average energy of the ionizing photons. Thus the higher the recombina- 
tion rate, the faster this effect will tend to push the temperature towards 

The average energy of the ionizing photons is given by the spectrum 
P{y) as {Euv) = h{v), where 

Here a is given by equation (|5|). Note that, in contrast to certain nebula 
calculations where all photons get absorbed sooner or later, the spectrum 
should be weighted by the photoionization cross section. This is because 
most photons never get absorbed, and all that is relevant is the energy 
distribution of those photons that do. Also note that P{v) is the energy 
distribution {W/ Hz)^ not the number distribution which is proportional to 
P{v)/v. 

The spectral parameters (Euv) and T* are given in Table I for some 
selected spectra. A power law spectrum P{iy) oc with a = 1 fits ob- 
served QSO spectra rather well in the vicinity of the Lyman limit (Ch- 
eney & Rowan-Robinson 1981; O'Brien et al. 1988), and is also consis- 
tent with the standard model for black hole accretion. A Planck spectrum 
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UV source 


Spectrum P{h>) 


(Euv) 




03 star 


T = 


50, OOOK Planck 


17.3 eV 


28,300K 


06 star 


T = 


40, OOOK Planck 


16.6 eV 


23,400K 


09 star 


T = 


30, OOOK Planck 


15.9 eV 


18,000K 


Pop. Ill star 


T = 


50, OOOK Vacca 


18.4 eV 


36,900K 


Black hole, QSO 


a = 


1 power law 


18.4 eV 


37,400K 


? 


a = 


2 power law 


17.2 eV 


27,800K 


? 


a = 


power law 


20.9 eV 


56,300K 


? 


T = 


100, OOOK Planck 


19.9 eV 


49,000K 



Table 1: Spectral parameters 



P(z^) oc z/^/ (^^^/^'^ — ij gives a decent prediction of T* for stars with sur- 
face temperatures below 30, OOOK. For very hot stars, more realistic spectra 
(Vacca 1993) fall off much slower above the Lyman limit, thus giving higher 
values of T* . As seen in Table |l], an extremely metal poor star of surface 
temperature 50, OOOK gives roughly the same T* as QSO radiation. The 
only stars that are likely to be relevant to early photoionization scenarios 
are extremely hot and short-lived ones, since the universe is less than a mil- 
lion years old at z = 1000, and fainter stars would be unable to inject enough 
energy in so short a time. Conceivably, less massive stars could play a the 
dominant role later on, thus lowering T* . However, since they radiate such 
a small fraction of their energy above the Lyman limit, very large numbers 
would be needed, which could be difficult to reconcile with the absence of 
observations of Population III stars today. 
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3.3 The thermal evolution 

At the low temperatures involved, the two dominant cooling effects^ are 
Compton drag against the microwave background photons and cooling due 
to the adiabatic expansion of the universe. Combining these effects, we 
obtain the following equation for the thermal evolution of the IGM: 



dT 

d{-z) 

where 



l + z 



-T + 



l + z 



T) + -Xrec{T){T* -T) 



(6) 



47r2 fkT. 



^comp 



45 



he J Home 



[l + z 



,-3 



0.00417/i-Vl + z) 



is (1 + z) 
r^o(l + z) 



^ times the Compton cooling rate per Hubble time and 



The factor of ^ in front of A,- 



is due to the fact that the 



photoelectrons end up sharing their energy with the protons. We have taken 
T^Q ~ 2.726K (Mather et al. 1994). Numerical solutions to this equation 
are shown in Figure ||, and the resulting y-parameters are given in Table ^. 

The temperature evolution separates into two distinct phases. In the first 
phase, which is almost instantaneous due to the high recombination rates at 
low temperatures, T rises very rapidly, up to a quasi-equilibrium tempera- 
ture slightly above the temperature of the microwave background photons. 
After this, in the second phase, T changes only slowly, and is approximately 
given by setting the expression in square brackets in equation (^) equal to 
zero. This quasi-equilibrium temperature is typically much lower than T*, 



^ Another cooling mechanism is coUisional excitation of atomic hydrogen followed by 
radiative de-excitation, which cools the IGM at a rate of (Dalgarno & McCray 1972) 



He 



7.5 X 10"^''e" 



"•8/T4„2 



n (1 — x)x erg cm s 



The ratio of this cooling rate to the Compton cooling rate is 



He 



exp [9.2 - 11.8/r4] (1 



(1 + Z)T4 



a quantity which is much smaller than unity for any reasonable parameter values when 
T < lO'^K. As will be seen, the temperatures at 2 ~ 1000 are typically a few thousand K, 
which with h^Qigm < 0.1 and a; > 0.9 renders coUisional excitation cooling more than nine 
orders of magnitude weaker than Compton cooling. Hence we can safely neglect coUisional 
excitations when computing the IGM temperature, the reason essentially being that the 
temperatures are so low that this process is suppressed by a huge Boltzman factor. 
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Model 


no 




h 




Zion 


1//0.000025 


Verdict 


QSO BDM I 


0.15 


0.15 


0.8 


37,400K 


1100 


6.30 


Ruled out 


QSO BDM II 


0.15 


0.15 


0.8 


37,400K 


200 


1.43 


Ruled out 


09 BDM 


0.15 


0.15 


0.8 


18,000K 


800 


2.91 


Ruled out 


QSO BDM III 


0.15 


0.04 


0.8 


37,400K 


1100 


0.67 


OK 


QSO CDM I 


1 


0.06 


0.5 


37,400K 


1100 


0.17 


OK 


QSO CDM II 


1 


0.03 


0.8 


37,400K 


1100 


0.16 


OK 



Table 2: Compton y-parameters for various scenarios 



since Compton cooling is so efficient at the high redshifts involved, and is 
given by 

AT = Te - P« - Te) OC -^g{Te)h^Qig^{T* - T), (7) 

independent of Qq, where g{Te) gS T"*^'^ encompasses the temperature de- 
pendence of Xrec- We typically have T <^T*. Using this, making the crude 
approximation of neglecting the temperature dependence of Xrec, and sub- 
stituting equation (^) into equation (|l]) indicates that 

Numerically selecting the best power-law fit, we find that this is indeed not 
too far off: the approximation 

y ^ omi2h^-'njgln,'/\T:f'{zion/noof' (8) 

is accurate to about 10% within the parameter range of cosmological interest. 
We have used equation (^) in Figure |^ by setting y = 2.5 x 10~^ and Zjon = 
1100. The shaded region of parameter space is thus ruled out by the COBE 
FIRAS experiment for fully ionized scenarios. 



4 Conclusions 

A reanalysis of the Compton y-distortion arising from early reionization 
shows that despite the radical sharpening of the FIRAS limit on y, one 
still cannot conclusively assert that the universe underwent a neutral phase. 
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Non-recombining scenarios where the ionization is thermal, caused by IGM 
temperatures exceeding lO^K, are clearly ruled out. Rather, the loophole 
is for the dominant ionization mechanism to be photoionization. We have 
shown that for spectra characteristic of both QSO radiation and massive 
metal-poor stars, the resulting IGM temperatures are so low that typical 
CDM models with no recombination can still survive the FIRAS test by a 
factor of six. This conclusion is valid if the flux of ionizing radiation is not 
so extreme that Compton heating becomes important. This is not difficult 
to arrange, as the cross section for Thomson scattering is some six orders of 
magnitude smaller than that for photoionization. 

For BDM models, the constraints are sharper. Non-recombining "clas- 
sical" BDM models with ^igm = ^^o ~ 0.15 are ruled out even with the 
extremely cautious reheating assumptions used in this paper, the earliest 
ionization redshift allowed being z ^ 130. Such models involving early 
non-linear seeds that on energetic grounds can very plausibly provide a pho- 
toionization source capable of reionizing the universe soon after the period of 
first recombination inevitably generate Compton distortions of order 10~^. 
These include texture as well as BDM models, both of which postulate, and 
indeed require, early reionization (z > 100) to avoid the generation of ex- 
cessive anisotropy in the cosmic microwave background on degree angular 
scales. 

Thus BDM models with reionization at z !=si 1000 can survive the y test 
only if most of the baryons form BDM when reionization occurs, and are 
thereby removed as a source of y-distortion, at least in the diffuse phase. 
This may be difficult to arrange at z > 100, since once the matter is reionizcd 
at this high a redshift, Compton drag is extremely effective in inhibiting any 
further gas collapse until z < 100. Since it takes only a small fraction of the 
baryons in the universe to provide a source of photons sufficient to maintain 
a fully ionized IGM even at z ^ 1000, we suspect that most of the baryons re- 
main diffuse until Compton drag eventually becomes ineffective. Moreover, 
the possibility that the IGM is only partially reionized at z ~ 1000 {e.g. 
GO), a situation which allows a lower value of the y-paramctcr, seems to us 
to be implausible as a delicate adjustment of ionization and recombination 
time-scales over a considerable range in z would be required. A comple- 
mentary argument that greatly restricts the parameter space allowable for 
fully ionized BDM models appeals to temperature fluctuations induced on 
the secondary last scattering surface, both by first order Doppler terms on 
degree scales and by second order terms on subarcminute scales (Hu et al. 
1994). Thus, BDM models would seem to be in some difficulty because of 
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the low limit on a possible y-distortion. 

The authors would like to thank W. Hu, A. Reisenegger and D. Scott 
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Figure 1: Thermal histories for various models. 
The temperature of the photoionized IGM is plotted for four of the cos- 
mological models and spectra of ionizing radiation listed in Table 2. The 
lowermost curve gives the temperature of the CMB photons. 



12 




1000 800 600 400 200 

1+z 



Figure 2: dy/dz for various models. 
The contribution to the y-parameter from different redshifts is plotted four 
of the cosmological models and spectra of ionizing radiation listed in Table 2. 
Thus for each model, the area under the curve is the predicted y-parameter. 
The area under the horizontal dashed line is 2.5 x 10~^, i.e. the COBE 
FIRAS limit. 
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Figure 3: Predicted and ruled out regions of parameter space. 
The hatched regions of parameter space are ruled out by the the COBE 
FIRAS limit y < 2.5 x 10"^ for Zion = 1100. J^o = 1 in the CDM plot 
and VIq = 0.15 in the BDM plot. The rectangular regions are the assumed 
parameter values for the CDM and BDM models, respectively. For CDM, 
the range 0.012 < h^l^VLigra < 0.024 is given by the nucleosynthesis con- 
straint 0.010 < h^^b < 0.015 and the assumption that 0.5 < h < 0.8. (If 
Qigm < ^b, the rectangle shifts to the left.) For the BDM models, h = 0.8 
and 0.03 < i^igm ^ ^o- The vertical range corresponds to feasible values of 
the spectral parameter T*. The upper limit corresponds to highly specula- 
tive star with surface temperature 100, OOOK and T* = 49, OOOK. The lower 
line corresponds to an 09 star. The dotted horizontal line corresponds to 
the spectrum expected from quasars/accreting black holes. 
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